Evolution and development of male-specific leg brushes in Drosophilidae.

The origin, diversification, and secondary loss of sexually dimorphic characters are common in animal evolution. In some cases, structurally and functionally similar traits have evolved independently in multiple lineages. Prominent examples of such traits include the male-specific grasping structures that develop on the front legs of many dipteran insects. In this report, we describe the evolution and development of one of these structures, the male-specific "sex brush." The sex brush is composed of densely packed, irregularly arranged modified bristles and is found in several distantly related lineages in the family Drosophilidae. Phylogenetic analysis using 250 genes from over 200 species provides modest support for a single origin of the sex brush followed by many secondary losses; however, independent origins of the sex brush cannot be ruled out completely. We show that sex brushes develop in very similar ways in all brush-bearing lineages. The dense packing of brush hairs is explained by the specification of bristle precursor cells at a near-maximum density permitted by the lateral inhibition mechanism, as well as by the reduced size of the surrounding epithelial cells. In contrast to the female and the ancestral male condition, where bristles are arranged in stereotypical, precisely spaced rows, cell migration does not contribute appreciably to the formation of the sex brush. The complex phylogenetic history of the sex brush can make it a valuable model for investigating coevolution of sex-specific morphology and mating behavior.

Diversity and function of arthropod endosymbiont toxins.

Bacterial endosymbionts induce dramatic phenotypes in their arthropod hosts, including cytoplasmic incompatibility, feminization, parthenogenesis, male killing, parasitoid defense, and pathogen blocking. The molecular mechanisms underlying these effects remain largely unknown but recent evidence suggests that protein toxins secreted by the endosymbionts play a role. Here, we describe the diversity and function of endosymbiont proteins with homology to known bacterial toxins. We focus on maternally transmitted endosymbionts belonging to the Wolbachia, Rickettsia, Arsenophonus, Hamiltonella, Spiroplasma, and Cardinium genera because of their ability to induce the above phenotypes. We identify at least 16 distinct toxin families with diverse enzymatic activities, including AMPylases, nucleases, proteases, and glycosyltransferases. Notably, several annotated toxins contain domains with homology to eukaryotic proteins, suggesting that arthropod endosymbionts mimic host biochemistry to manipulate host physiology, similar to bacterial pathogens.

Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola.

Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.

Co-evolving wing spots and mating displays are genetically separable traits in Drosophila.

The evolution of sexual traits often involves correlated changes in morphology and behavior. For example, in Drosophila, divergent mating displays are often accompanied by divergent pigment patterns. To better understand how such traits co-evolve, we investigated the genetic basis of correlated divergence in wing pigmentation and mating display between the sibling species Drosophila elegans and Drosophila gunungcola. Drosophila elegans males have an area of black pigment on their wings known as a wing spot and appear to display this spot to females by extending their wings laterally during courtship. By contrast, D. gunungcola lost both of these traits. Using Multiplexed Shotgun Genotyping (MSG), we identified a ∼440 kb region on the X chromosome that behaves like a genetic switch controlling the presence or absence of male-specific wing spots. This region includes the candidate gene optomotor-blind (omb), which plays a critical role in patterning the Drosophila wing. The genetic basis of divergent wing display is more complex, with at least two loci on the X chromosome and two loci on autosomes contributing to its evolution. Introgressing the X-linked region affecting wing spot development from D. gunungcola into D. elegans reduced pigmentation in the wing spots but did not affect the wing display, indicating that these are genetically separable traits. Consistent with this observation, broader sampling of wild D. gunungcola populations confirmed that the wing spot and wing display are evolving independently: some D. gunungcola males performed wing displays similar to D. elegans despite lacking wing spots. These data suggest that correlated selection pressures rather than physical linkage or pleiotropy are responsible for the coevolution of these morphological and behavioral traits. They also suggest that the change in morphology evolved prior to the change in behavior.

The yellow gene influences Drosophila male mating success through sex comb melanization.

Drosophila melanogaster males perform a series of courtship behaviors that, when successful, result in copulation with a female. For over a century, mutations in the yellow gene, named for its effects on pigmentation, have been known to reduce male mating success. Prior work has suggested that yellow influences mating behavior through effects on wing extension, song, and/or courtship vigor. Here, we rule out these explanations, as well as effects on the nervous system more generally, and find instead that the effects of yellow on male mating success are mediated by its effects on pigmentation of male-specific leg structures called sex combs. Loss of yellow expression in these modified bristles reduces their melanization, which changes their structure and causes difficulty grasping females prior to copulation. These data illustrate why the mechanical properties of anatomy, not just neural circuitry, must be considered to fully understand the development and evolution of behavior.

Pleiotropic Effects of ebony and tan on Pigmentation and Cuticular Hydrocarbon Composition in Drosophila melanogaster.

Pleiotropic genes are genes that affect more than one trait. For example, many genes required for pigmentation in the fruit fly Drosophila melanogaster also affect traits such as circadian rhythms, vision, and mating behavior. Here, we present evidence that two pigmentation genes, ebony and tan, which encode enzymes catalyzing reciprocal reactions in the melanin biosynthesis pathway, also affect cuticular hydrocarbon (CHC) composition in D. melanogaster females. More specifically, we report that ebony loss-of-function mutants have a CHC profile that is biased toward long (>25C) chain CHCs, whereas tan loss-of-function mutants have a CHC profile that is biased toward short (<25C) chain CHCs. Moreover, pharmacological inhibition of dopamine synthesis, a key step in the melanin synthesis pathway, reversed the changes in CHC composition seen in ebony mutants, making the CHC profiles similar to those seen in tan mutants. These observations suggest that genetic variation affecting ebony and/or tan activity might cause correlated changes in pigmentation and CHC composition in natural populations. We tested this possibility using the Drosophila Genetic Reference Panel (DGRP) and found that CHC composition covaried with pigmentation as well as levels of ebony and tan expression in newly eclosed adults in a manner consistent with the ebony and tan mutant phenotypes. These data suggest that the pleiotropic effects of ebony and tan might contribute to covariation of pigmentation and CHC profiles in Drosophila.

Author Correction: Genotypic variability-based genome-wide association study identifies non-additive loci HLA-C and IL12B for psoriasis.

This article was originally published under a CC BY-NC-SA License, but has now been made available under a CC BY 4.0 License.

Genome-wide association study of response to tumour necrosis factor inhibitor therapy in rheumatoid arthritis.

Rheumatoid arthritis (RA) is characterised by chronic synovial joint inflammation. Treatment has been revolutionised by tumour necrosis factor alpha inhibitors (TNFi) but each available drug shows a significant non-response rate. We conducted a genome-wide association study of 1752 UK RA TNFi-treated patients to identify predictors of change in the Disease Activity Score 28 (DAS28) and subcomponents over 3-6 months. The rs7195994 variant at the FTO gene locus was associated with infliximab response when looking at a change in the swollen joint count (SJC28) subcomponent (p = 9.74 × 10-9). Capture Hi-C data show chromatin interactions in GM12878 cells between rs2540767, in high linkage disequilibrium with rs7195994 (R2 = 0.9) and IRX3, a neighbouring gene of FTO. IRX3 encodes a transcription factor involved in adipocyte remodelling and is regarded as the obesity gene at the FTO locus. Importantly, the rs7195994 association remained significantly associated following adjustment for BMI. In addition, using capture Hi-C data we showed interactions between TNFi-response associated variants and 16 RA susceptibility variants.

Genome-wide association study of response to methotrexate in early rheumatoid arthritis patients.

Methotrexate (MTX) monotherapy is a common first treatment for rheumatoid arthritis (RA), but many patients do not respond adequately. In order to identify genetic predictors of response, we have combined data from two consortia to carry out a genome-wide study of response to MTX in 1424 early RA patients of European ancestry. Clinical endpoints were change from baseline to 6 months after starting treatment in swollen 28-joint count, tender 28-joint count, C-reactive protein and the overall 3-component disease activity score (DAS28). No single nucleotide polymorphism (SNP) reached genome-wide statistical significance for any outcome measure. The strongest evidence for association was with rs168201 in NRG3 (p = 10-7 for change in DAS28). Some support was also seen for association with ZMIZ1, previously highlighted in a study of response to MTX in juvenile idiopathic arthritis. Follow-up in two smaller cohorts of 429 and 177 RA patients did not support these findings, although these cohorts were more heterogeneous.

Evaluation of metrics for benchmarking antimicrobial use in the UK dairy industry.

The issue of antimicrobial resistance is of global concern across human and animal health. In 2016, the UK government committed to new targets for reducing antimicrobial use (AMU) in livestock. Although a number of metrics for quantifying AMU are defined in the literature, all give slightly different interpretations. This paper evaluates a selection of metrics for AMU in the dairy industry: total mg, total mg/kg, daily dose and daily course metrics. Although the focus is on their application to the dairy industry, the metrics and issues discussed are relevant across livestock sectors. In order to be used widely, a metric should be understandable and relevant to the veterinarians and farmers who are prescribing and using antimicrobials. This means that clear methods, assumptions (and possible biases), standardised values and exceptions should be published for all metrics. Particularly relevant are assumptions around the number and weight of cattle at risk of treatment and definitions of dose rates and course lengths; incorrect assumptions can mean metrics over-represent or under-represent AMU. The authors recommend that the UK dairy industry work towards the UK-specific metrics using the UK-specific medicine dose and course regimens as well as cattle weights in order to monitor trends nationally.

Genotypic variability-based genome-wide association study identifies non-additive loci HLA-C and IL12B for psoriasis.

Genome-wide association studies (GWASs) have identified a number of loci for psoriasis but largely ignored non-additive effects. We report a genotypic variability-based GWAS (vGWAS) that can prioritize non-additive loci without requiring prior knowledge of interaction types or interacting factors in two steps, using a mixed model to partition dichotomous phenotypes into an additive component and non-additive environmental residuals on the liability scale and then the Levene's (Brown-Forsythe) test to assess equality of the residual variances across genotype groups genome widely. The vGWAS identified two genome-wide significant (P < 5.0e-08) non-additive loci HLA-C and IL12B that were also genome-wide significant in an accompanying GWAS in the discovery cohort. Both loci were statistically replicated in vGWAS of an independent cohort with a small sample size. HLA-C and IL12B were reported in moderate gene-gene and/or gene-environment interactions in several occasions. We found a moderate interaction with age-of-onset of psoriasis, which was replicated indirectly. The vGWAS also revealed five suggestive loci (P < 6.76e-05) including FUT2 that was associated with psoriasis with environmental aspects triggered by virus infection and/or metabolic factors. Replication and functional investigation are needed to validate the suggestive vGWAS loci.

CD4+ and B Lymphocyte Expression Quantitative Traits at Rheumatoid Arthritis Risk Loci in Patients With Untreated Early Arthritis: Implications for Causal Gene Identification.

OBJECTIVE: Rheumatoid arthritis (RA) is a genetically complex disease of immune dysregulation. This study sought to gain further insight into the genetic risk mechanisms of RA by conducting an expression quantitative trait locus (eQTL) analysis of confirmed genetic risk loci in CD4+ T cells and B cells from carefully phenotyped patients with early arthritis who were naive to therapeutic immunomodulation. METHODS: RNA and DNA were isolated from purified B and/or CD4+ T cells obtained from the peripheral blood of 344 patients with early arthritis. Genotyping and global gene expression measurements were carried out using Illumina BeadChip microarrays. Variants in linkage disequilibrium (LD) with non-HLA RA single-nucleotide polymorphisms (defined as r2 ≥ 0.8) were analyzed, seeking evidence of cis- or trans-eQTLs according to whether the associated probes were or were not within 4 Mb of these LD blocks. RESULTS: Genes subject to cis-eQTL effects that were common to both CD4+ and B lymphocytes at RA risk loci were FADS1, FADS2, BLK, FCRL3, ORMDL3, PPIL3, and GSDMB. In contrast, those acting on METTL21B, JAZF1, IKZF3, and PADI4 were unique to CD4+ lymphocytes, with the latter candidate risk gene being identified for the first time in this cell subset. B lymphocyte-specific eQTLs for SYNGR1 and CD83 were also found. At the 8p23 BLK-FAM167A locus, adjacent genes were subject to eQTLs whose activity differed markedly between cell types; in particular, the FAM167A effect displayed striking B lymphocyte specificity. No trans-eQTLs approached experiment-wide significance, and linear modeling did not identify a significant influence of biologic covariates on cis-eQTL effect sizes. CONCLUSION: These findings further refine the understanding of candidate causal genes in RA pathogenesis, thus providing an important platform from which downstream functional studies, directed toward particular cell types, may be prioritized.

Analysis of SNP-SNP interactions and bone quantitative ultrasound parameter in early adulthood.

BACKGROUND: Osteoporosis individual susceptibility is determined by the interaction of multiple genetic variants and environmental factors. The aim of this study was to conduct SNP-SNP interaction analyses in candidate genes influencing heel quantitative ultrasound (QUS) parameter in early adulthood to identify novel insights into the mechanism of disease. METHODS: The study population included 575 healthy subjects (mean age 20.41; SD 2.36). To assess bone mass QUS was performed to determine Broadband ultrasound attenuation (BUA, dB/MHz). A total of 32 SNPs mapping to loci that have been characterized as genetic markers for QUS and/or BMD parameters were selected as genetic markers in this study. The association of all possible SNP pairs with QUS was assessed by linear regression and a SNP-SNP interaction was defined as a significant departure from additive effects. RESULTS: The pairwise SNP-SNP analysis showed multiple interactions. The interaction comprising SNPs rs9340799 and rs3736228 that map in the ESR1 and LRP5 genes respectively, revealed the lowest p value after adjusting for confounding factors (p-value = 0.001, ß (95% CI) = 14.289 (5.548, 23.029). In addition, our model reported others such as TMEM135-WNT16 (p = 0.007, ß(95%CI) = 9.101 (2.498, 15.704), ESR1-DKK1 (p = 0.012, ß(95%CI) = 13.641 (2.959, 24.322) or OPG-LRP5 (p = 0.012, ß(95%CI) = 8.724 (1.936, 15.512). However, none of the detected interactions remain significant considering the Bonferroni significance threshold for multiple testing (p<0.0001). CONCLUSION: Our analysis of SNP-SNP interaction in candidate genes of QUS in Caucasian young adults reveal several interactions, especially between ESR1 and LRP5 genes, that did not reach statistical significance. Although our results do not support a relevant genetic contribution of SNP-SNP epistatic interactions to QUS in young adults, further studies in larger independent populations would be necessary to support these preliminary findings.

Identifying Causal Genes at the Multiple Sclerosis Associated Region 6q23 Using Capture Hi-C.

BACKGROUND: The chromosomal region 6q23 has been found to be associated with multiple sclerosis (MS) predisposition through genome wide association studies (GWAS). There are four independent single nucleotide polymorphisms (SNPs) associated with MS in this region, which spans around 2.5 Mb. Most GWAS variants associated with complex traits, including these four MS associated SNPs, are non-coding and their function is currently unknown. However, GWAS variants have been found to be enriched in enhancers and there is evidence that they may be involved in transcriptional regulation of their distant target genes through long range chromatin looping. AIM: The aim of this work is to identify causal disease genes in the 6q23 locus by studying long range chromatin interactions, using the recently developed Capture Hi-C method in human T and B-cell lines. Interactions involving four independent associations unique to MS, tagged by rs11154801, rs17066096, rs7769192 and rs67297943 were analysed using Capture Hi-C Analysis of Genomic Organisation (CHiCAGO). RESULTS: We found that the pattern of chromatin looping interactions in the MS 6q23 associated region is complex. Interactions cluster in two regions, the first involving the rs11154801 region and a second containing the rs17066096, rs7769192 and rs67297943 SNPs. Firstly, SNPs located within the AHI1 gene, tagged by rs11154801, are correlated with expression of AHI1 and interact with its promoter. These SNPs also interact with other potential candidate genes such as SGK1 and BCLAF1. Secondly, the rs17066096, rs7769192 and rs67297943 SNPs interact with each other and with immune-related genes such as IL20RA, IL22RA2, IFNGR1 and TNFAIP3. Finally, the above-mentioned regions interact with each other and therefore, may co-regulate these target genes. CONCLUSION: These results suggest that the four 6q23 variants, independently associated with MS, are involved in the regulation of several genes, including immune genes. These findings could help understand mechanisms of disease and suggest potential novel therapeutic targets.

Capture Hi-C identifies a novel causal gene, IL20RA, in the pan-autoimmune genetic susceptibility region 6q23.

BACKGROUND: The identification of causal genes from genome-wide association studies (GWAS) is the next important step for the translation of genetic findings into biologically meaningful mechanisms of disease and potential therapeutic targets. Using novel chromatin interaction detection techniques and allele specific assays in T and B cell lines, we provide compelling evidence that redefines causal genes at the 6q23 locus, one of the most important loci that confers autoimmunity risk. RESULTS: Although the function of disease-associated non-coding single nucleotide polymorphisms (SNPs) at 6q23 is unknown, the association is generally assigned to TNFAIP3, the closest gene. However, the DNA fragment containing the associated SNPs interacts through chromatin looping not only with TNFAIP3, but also with IL20RA, located 680 kb upstream. The risk allele of the most likely causal SNP, rs6927172, is correlated with both a higher frequency of interactions and increased expression of IL20RA, along with a stronger binding of both the NFκB transcription factor and chromatin marks characteristic of active enhancers in T-cells. CONCLUSIONS: Our results highlight the importance of gene assignment for translating GWAS findings into biologically meaningful mechanisms of disease and potential therapeutic targets; indeed, monoclonal antibody therapy targeting IL-20 is effective in the treatment of rheumatoid arthritis and psoriasis, both with strong GWAS associations to this region.

Previously reported PDE3A-SLCO1C1 genetic variant does not correlate with anti-TNF response in a large UK rheumatoid arthritis cohort.

AIM: A genetic variant has recently reached genome-wide significance for association with TNF-inhibitor response in rheumatoid arthritis patients. Here we undertake a replication study in a UK Caucasian population to test for association with TNF-inhibitor response. MATERIALS & METHODS: The genetic variant, rs3794271, located within the PDE3A-SLCO1C1 locus was analyzed for correlation with treatment response using both the EULAR classification criteria and absolute change in (Δ)DAS28 scores as outcome measures. RESULTS: Genotype data were available from 1750 TNF-inhibitor treated individuals. However, no evidence for association was observed (EULAR: p = 0.91 and ΔDAS28: p = 0.93). Furthermore, no significant associations were observed upon stratification by the anti-TNF received (p > 0.05). CONCLUSION: In the largest replication cohort conducted to date, no evidence for association was observed.

Major histocompatibility complex harbors widespread genotypic variability of non-additive risk of rheumatoid arthritis including epistasis.

Genotypic variability based genome-wide association studies (vGWASs) can identify potentially interacting loci without prior knowledge of the interacting factors. We report a two-stage approach to make vGWAS applicable to diseases: firstly using a mixed model approach to partition dichotomous phenotypes into additive risk and non-additive environmental residuals on the liability scale and secondly using the Levene's (Brown-Forsythe) test to assess equality of the residual variances across genotype groups per marker. We found widespread significant (P < 2.5e-05) vGWAS signals within the major histocompatibility complex (MHC) across all three study cohorts of rheumatoid arthritis. We further identified 10 epistatic interactions between the vGWAS signals independent of the MHC additive effects, each with a weak effect but jointly explained 1.9% of phenotypic variance. PTPN22 was also identified in the discovery cohort but replicated in only one independent cohort. Combining the three cohorts boosted power of vGWAS and additionally identified TYK2 and ANKRD55. Both PTPN22 and TYK2 had evidence of interactions reported elsewhere. We conclude that vGWAS can help discover interacting loci for complex diseases but require large samples to find additional signals.

Replication of Associations of Genetic Loci Outside the HLA Region With Susceptibility to Anti-Cyclic Citrullinated Peptide-Negative Rheumatoid Arthritis.

OBJECTIVE: Genetic polymorphisms within the HLA region explain only a modest proportion of anti-cyclic citrullinated peptide (anti-CCP)-negative rheumatoid arthritis (RA) heritability. However, few non-HLA markers have been identified so far. This study was undertaken to replicate the associations of anti-CCP-negative RA with non-HLA genetic polymorphisms demonstrated in a previous study. METHODS: The Rheumatoid Arthritis Consortium International densely genotyped 186 autoimmune-related regions in 3,339 anti-CCP-negative RA patients and 15,870 controls across 6 different populations using the Illumina ImmunoChip array. We performed a case-control replication study of the anti-CCP-negative markers with the strongest associations in that discovery study, in an independent cohort of anti-CCP-negative UK RA patients. Individuals from the arcOGEN Consortium and Wellcome Trust Case Control Consortium were used as controls. Genotyping in cases was performed using Sequenom MassArray technology. Genome-wide data from controls were imputed using the 1000 Genomes Phase I integrated variant call set release version 3 as a reference panel. RESULTS: After genotyping and imputation quality control procedures, data were available for 15 non-HLA single-nucleotide polymorphisms in 1,024 cases and 6,348 controls. We confirmed the known markers ANKRD55 (meta-analysis odds ratio [OR] 0.80; P = 2.8 × 10(-13) ) and BLK (OR 1.13; P = 7.0 × 10(-6) ) and identified new and specific markers of anti-CCP-negative RA (prolactin [PRL] [OR 1.13; P = 2.1 × 10(-6) ] and NFIA [OR 0.85; P = 2.5 × 10(-6) ]). Neither of these loci is associated with other common, complex autoimmune diseases. CONCLUSION: Anti-CCP-negative RA and anti-CCP-positive RA are genetically different disease subsets that only partially share susceptibility factors. Genetic polymorphisms located near the PRL and NFIA genes represent examples of genetic susceptibility factors specific for anti-CCP-negative RA.

A role for HLA-DRB1*1101 and DRB1*0801 in cognitive ability and its decline with age.

Cognitive abilities (memory, processing speed, vocabulary, and fluid intelligence) are correlated with educational attainment and occupational status, as well as physical and mental health. The variation in cognitive abilities observed within a population has a substantial genetic contribution (heritability ∼50%) and yet the identification of genetic polymorphisms from both genome-wide association and candidate studies have to date only uncovered a limited number of genetic variants that exert small genetic effects. Here we impute human leukocyte antigens (HLA) using existing genome-wide association data from 1,559 non-pathological elderly volunteers who have been followed for changes in cognitive functioning between a 12- and 18-year period. Specifically, we investigate DRB1*05 (*11/*12) and DRB1*01, which have previously been associated with cognitive ability. We also analyze DRB1*0801, which shares close sequence homology with DRB1*1101. Together with DRB1*1101, DRB1*0801 has been associated with several diseases including multiple sclerosis and primary biliary cirrhosis, which themselves are associated with cognitive impairment. We observed that both DRB1*0801 and DRB1*1101 were significantly associated with vocabulary ability (cross-sectional and longitudinal scores) and that the effects were in opposite directions with DRB1*0801 associated with lower score and faster decline. This opposing affect is similar to that reported by other groups in systemic lupus erythematosus, type 1 diabetes, and primary biliary cirrhosis. DRB1*0801 was also significantly associated with reduced memory ability. We observed no associations between cognitive abilities and DRB1*01 or DRB1*12.